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Abstract 

We address the following problem: Given a complete fc-partite geometric graph K whose 
vertex set is a set of n points in R d , compute a spanner of K that has a "small" stretch factor 
and "few" edges. We present two algorithms for this problem. The first algorithm computes 
a (5 + e)-spanner of K with 0(n) edges in 0(n log n) time. The second algorithm computes a 
(3 + e)-spanner of K with O(nlogn) edges in 0(n log n) time. The latter result is optimal: We 
show that for any 2 < k < n — 8(yn log n), spanners with 0(n log n) edges and stretch factor 
less than 3 do not exist for all complete fc-partite geometric graphs. 

1 Introduction 

Let S be a set of n points in A geometric graph with vertex set S is an undirected graph H 
whose edges are line segments pq that are weighted by the Euclidean distance \pq\ between p and 
q. For any two points p and q in S, we denote by Sn(p,q) the length of a shortest path in H 
between p and q. For a real number t > 1, a subgraph G of H is said to be a t- spanner of H, if 
3g(p,q) < t • 5h(p,q) for all points p and q in S. The smallest i for which this property holds is 
called the stretch factor of G. Thus, a subgraph G of H with stretch factor t approximates the Q) 
pairwise shortest-path lengths in if within a factor of i. If H is the complete geometric graph with 
vertex set S, then G is also called a i-spanner of the point set S. 

Most of the work on constructing spanners has been done for the case when H is the complete 
graph. It is well known that for any set S of n points in W 1 and for any real constant e > 0, there 
exists a (1 + e)-spanner of S containing 0{n) edges. Moreover, such spanners can be computed in 
0{n log n) time; see Salowe [H] and Vaidya jS]. For a detailed overview of results on spanners for 
point sets, see the book by Narasimhan and Smid [6]. 

For spanners of arbitrary geometric graphs, much less is known. Althofer et al. [1] have shown 
that for any t > 1, every weighted graph H with n vertices contains a subgraph with 0(n 1+2 ^* -1 ^) 
edges, which is a t-spanner of H. Observe that this result holds for any weighted graph; in particular, 
it is valid for any geometric graph. For geometric graphs, a lower bound was given by Gudmundsson 
and Smid [5]: They proved that for every real number t with 1 < t < | log n, there exists a geometric 
graph H with n vertices, such that every t-spanner of H contains 0(n 1+1// *) edges. Thus, if we are 
looking for spanners with 0{n) edges of arbitrary geometric graphs, then the best stretch factor 
we can obtain is O(logn). 
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In this paper, we consider the case when the input graph is a complete /c-partite geometric 
graph. Let S be a set of n points in M. d , and let S be partitioned into subsets Ci, C2, . . . , Cfc. Let 
Kcx...c k denote the complete k-partite graph on S. This graph has 5 as its vertex set and two points 
p and q are connected by an edge (of length \pq\) if and only if p and q are in different subsets of 
the partition. The problem we address is formally defined as follows: 

Problem 1.1 Let k > 2 be an integer, let S be a set of n points in M. d , and let S be partitioned 
into k subsets C%,C2, ■ ■ ■ ,Cj~. Compute a t-spanner of the complete k-partite graph Kc x ...c k that 
has a "small" number of edges and whose stretch factor t is "small". 

The main contribution of this paper is to present an algorithm that computes such a t-spanner 
with 0{n) edges in O(nlogn) time, where t = 5 + e for any constant e > 0. We also show that if 
one is willing to use 0{n log n) edges, then our algorithm adapts easily to reach a stretch factor of 
t = 3 + e. Finally, we show that the latter result is optimal: For any k with 2 < k < n — 0(^/nlogn), 
spanners with 0{n log n) edges and stretch factor less than 3 do not exist for all complete fc-partite 
geometric graphs. 

We remark that in a recent paper, Bose et al. (2] considered the problem of constructing spanners 
of point sets that have 0{n) edges and whose chromatic number is a most k. This problem is 
different from ours: Bose et al. compute a spanner of the complete graph and their algorithm can 
choose a "good" fc-partition of the vertices. In our problem, the fc-partition is given and we want 
to compute a spanner of the complete /c-partite graph. 

Possible applications of our algorithm are in wireless networks having the property that com- 
municating nodes are partitioned into sets such that two nodes can communicate if and only if 
they do not belong to the same set. This would be the case, for example, when Time Division 
Multiplexing (TDMA) is used. Since the wireless medium prohibits simultaneous transmission and 
reception at one node, two nodes communicating during the same time slots cannot communicate 
with each other; see Raman and Chebrolu [7|. 

The rest of this paper is organized as follows. In Section [2] we recall properties of the Well- 
Separated Pair Decomposition (WSPD) that we use in our algorithm. In Section |3j we provide 
an algorithm that solves the problem of constructing a spanner of the complete /c-partite graph. 
In Section [4] we show that the spanner constructed by this algorithm has 0(n) edges and that 
its stretch factor is bounded from above by a constant that depends only on the dimension d. In 
Section [5] we show how a simple modification to our algorithm improves the stretch factor to 5 + e 
while still having 0(n) edges. In Section [6] we show how to achieve a stretch factor of 3 + e using 
0(n log n) edges. We also prove that the latter result is optimal. We conclude in Section [7] 

2 The Well- Separated Pair Decomposition 

In this section, we recall crucial properties of the Well-Separated Pair Decomposition (WSPD) of 
Callahan and Kosaraju jl] that we use for our construction. The reader who is familiar with the 
WSPD may go directly to Section [3] Our presentation follows the one in Narasimhan and Smid [6]. 
Intuitively, a WSPD is a partition of the edges of a complete geometric graph such that all edges 
that are grouped together are approximately equal. To give a formal definition of the WSPD, we 
first need to define what it means for two sets to be well-separated. 
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Definition 2.1 Let S be a set of points in W 1 . The bounding box @(S) of S is the smallest 
axes-parallel hyperrectangle that contains S. 

Definition 2.2 Let X and Y be two sets of points in M d and let s > be a real number. We say 
that X and Y are well-separated with respect to s if there exists two balls B\ and L>2 such that 

1. B\ and B<i have the same radius, say p, 

2. B\ contains the bounding box of X, 

3. i?2 contains the bounding box ofY, and 

4- the distance mm{\xy\ : x G B\ n M. d , y G B2 n M d } between B\ and B2 is at least sp. 

Definition 2.3 Let S be a set of points in It d and let s > be a real number. A well-separated 
pair decomposition (WSPD) of S with separation constant s is a set of unordered pairs of subsets 
of S that are well-separated with respect to s, such that for any two distinct points p,q G S there is 
a unique pair {X, Y} in the WSPD such that p G X and q £Y. 

Lemma 2.4 (Lemma 9.1.2 in |6j) Let s > be a real number and let X and Y be two point sets 
that are well-separated with respect to s. 

1. Ifp,p',p" G X and q £Y, then \p'p"\ < (2/s)\pq\. 

2. Ifp,p' G X and q,q' G Y, then \p'q'\ < (1 + 4/s)\pq\. 

The first part of this lemma states that distances within one set are very small compared to 
distances between pairs of points having one endpoint in each set. The second part states that all 
pairs of points having one endpoint in each set have approximately the same distance. 

Callahan and Kosaraju [3] have shown how to construct a t-spanner of S from a WSPD: All 
one has to do is pick from each pair {X, Y} an arbitrary edge (p, q) with and q £Y. Using 

induction on the rank of the length of the edges in the complete graph K$, it can be shown that, 
when s > 4, this process leads to a ((s + 4)/(s — 4))-spanner. Thus, by choosing s to be a sufficiently 
large constant, the stretch factor can be made arbitrarily close to 1. 

In order to compute a spanner of S that has a linear number of edges, one needs a WSPD that 
has a linear number of pairs. Callahan and Kosaraju [4 showed that a WSPD with a linear number 
of pairs always exists and can be computed in time 0{n log n). Their algorithm uses a split-tree. 

Definition 2.5 Let S be a non-empty set of points in W*. The split-tree of S is defined as follows: 
if S contains only one point, then the split-tree is a single node that stores that point. Otherwise, 
the split-tree has a root that stores the bounding box (3(S) of S, as well as an arbitrary point of S 
called the representative of S and denoted by rep(S). Split f3(S) into two hyperrectangles by cutting 
its longest interval into two equal parts, and let Si and S2 be the subsets of S contained in the 
two hyperrectangles. The root of the split-tree of S has two sub-trees, which are recursively defined 
split-trees of S\ and S^- 

The precise way Callahan and Kosaraju used the split-tree to compute a WSPD with a linear 
number of pairs is of no importance to us. The only important aspect we need to retain is that 
each pair is uniquely determined by a pair of nodes in the tree. More precisely, for each pair {X, Y} 
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in the WSPD that is output by their algorithm, there are unique internal nodes u and v in the 
split-tree such that the sets S u and S v of points stored at the leaves of the subtrees rooted at u 
and v are precisely X and Y. Since there is such a unique correspondence, we will denote pairs in 
the WSPD by {S u , S v }, meaning that u and v are the nodes corresponding to the sets X = S u and 
Y = S v . Also, although the WSPD of a point set is not unique, when we talk about the WSPD, 
we mean the WSPD that is computed by the algorithm of Callahan and Kosaraju. 

Before we present our algorithm, we give the statement of the following lemmas that we use 
to analyze our algorithm in Section |4j If R is an axes-parallel hyperrectangle in then we use 
L max .(R) to denote the length of a longest side of R. 

Lemma 2.6 (Lemma 9.5.3 in |6j) Let u be a node in the split-tree and let u' be a node in the 
subtree of u such that the path between them contains at least d edges. Then 

LraMSu')) < \ ■ L max (P(S u )). 

Lemma 2.7 (Lemma 11.3.1 in [6j) Let {S U ,S V } be a pair in the WSPD, let t be the distance 
between the centers of f3(S u ) and /3(S V ), and let n(u) be the parent of u in the split-tree. Then 

21 

Vd(s + 4) 

3 A First Algorithm 

We now show how the WSPD can be used to address the problem of computing a spanner of a 
complete fc-partite graph. In this section, we introduce an algorithm that outputs a graph with 
constant stretch factor and 0(n) edges. The analysis of this algorithm is presented in Section |4j 
In Section [5] we show how this algorithm can be improved to achieve a stretch factor of 5 + e. 

The input set S C M. d is the disjoint union of k sets C\, C2, . . . , Cp,. We say that the elements 
of C c have "color" c. The graph K = Kc 1 ...c k is the complete fc-partite geometric graph. 

Definition 3.1 Let T be the split-tree of S that is used to compute the WSPD of S. 

1. For any node u in T, we denote by S u the set of all points in the subtree rooted at u. 

2. We define MWSPD to be the subset of the WSPD obtained by removing all pairs {S u , S v } for 
which all points of S u U S v have the same color. 

3. A node u inT is called multichromatic if there exist points p and q in S u and a node v in T, 
such that p and q have different colors and {S u , S v } is in the MWSPD. 

4- A node u in T is called a c-node if all points of S u have color c and there exists a node v in 
T such that {S u , S v } is in the MWSPD. 

5. A c-node u in T is called a c-root if it does not have a proper ancestor that is a c-node in T. 

6. A c-node u inT is called a c-leaf if it does not have another c-node in its subtree. 

7. A c-node v! in T is called a c-child of a c-node u inT if u' is in the subtree rooted at u and 
there is no c-node on the path strictly between u and v! . 
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All points are of the same color. 




At least one set only contains points of one color 
but not all points are of the same color. 

Figure 1: The three cases of Algorithm [T] 

8. For each color c and for each c-node u in T, rep{S u ) denotes a fixed arbitrary point in S u . 

9. For each multichromatic node u in T, rep{S u ) and rep'{S u ) denote two fixed arbitrary points 
in S u that have different colors. 

10. The distance between two sets S v and S w , denoted by dist(S v , S w ), is defined to be the distance 
between the centers of their bounding boxes. 

11. Let u be a c-node in T. Consider all pairs {S v ,Sw\ in the MWSPD, where v is a c-node on 
the path in T from u to the root (this path includes u). Let {S v , S w } be such a pair for which 
dist(S v , S w ) is minimum. We define cl(S u ) to be the set S w . 

Algorithm [T] computes a spanner of a complete fc-partite geometric graph K = Kc x ...c k - The 
intuition behind this algorithm is the following. First, the algorithm computes the WSPD. Then, 
it considers each pair {S u , S v } of the WSPD, and decides whether or not to add an edge between 
S u and S v . The outcome of this decision is based on the following three cases, which are illustrated 
in Figures [T] and [2j 

Case 1: All points of S u U S v are of the same color. In this case, there is no edge of K to 
approximate, so the algorithm ignores this pair. 

Case 2: Both S u and S v are multichromatic. In this case, the algorithm adds one edge between 
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Figure 2: Handling a c-node. 



S u and S v to the spanner; see lines 28 - 29 Observe that the two vertices of this edge do not have 
the same color. This edge will allow us to approximate each edge (p, q) of K, where p S S u , q £ S v , 
and p and q have different colors. 

Case 3: All points in S u are of the same color c. In this case, an edge is added between iep(S u ) and 
one of the two representatives of S v whose color is not c; see lines 17-[T8} In order to approximate 
each edge of K having one vertex (of color c) in S u and the other vertex (of a different color) in S v , 
more edges have to be added. This is done in such a way that our final graph contains a "short" 
path between every point p of S u and the representative iep(S u ) of S u . Observe that this path 
must contain points whose color is not equal to c; thus, these points are not in S u . One way to 
achieve this is to add an edge between each point of S u and one of the two representatives of cl(S u ) 
whose color is not c; we call this construction a star. However, since the subtree rooted at u may 
contain other c-nodes, many edges may be added for each point in S u , which could possibly lead to 
a quadratic number of edges in the final graph. To guarantee that the algorithm does not add too 
many edges, it introduces a star only if u is a c-leaf; see lines [8f|ll[ If u is a c-node, the algorithm 
only adds one edge between rep(S u ) and a representatives of cl(S u ) whose color is not c; see lines 
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whose color is not c; see lines 



4 Analysis of Algorithm [T] 

Lemma 4.1 The graph G computed by Algorithm^has 0(\S\) edges. 
Proof: For each color c and for each c-leaf u', the algorithm adds \S u i\ edges to G in lines [9 10 



Since the sets S u i, where u' ranges over all c- leaves and c ranges over all colors, are pairwise disjoint, 
the total number of edges that are added in lines gjlO] is 0(\S\). 

The total number of edges that are added in lines |17fjl8| and [28] - p9| is at most the number of 
pairs in the MWSPD. Since the WSPD contains 0(|5|) pairs (see the same upper bound holds 



for the number of edges added in lines 17-18 and|28||29 



The total number of edges that are added in lines 14 - 15 and 21-22 is at most twice the number 
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Algorithm 1: Computing a sparse subgraph of Kc 1 ...c k whose stretch factor is bounded by 
a constant. 

Input: A set S of points in R d , which is partitioned into k subsets C±, . . . , C^. 

Output: A spanner G = (S, E) of the complete fe-partite graph Kc 1 ...c k - 

1 compute the split-tree T of S; 

2 using T, compute the WSPD with respect to a separation constant s > 0; 

3 using the WSPD, compute the MWSPD; 

5 for each color c in {1, 2, ... , k} do 

6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 



for each c-root u in T do 

for each c-leaf u' in the subtree of u do 
for each p £ S u / do 

if rep(cl(5 u ')) does not have color c then add (p,Tep(cl(S u >))) to E; 
else add (p, rep'(cl(<!v))) to E 1 ; 
end 
end 

for each c-node u' that is in the subtree of u (including u) do 

if rep(cl(<S u /)) does not have color c then add (rep(S' u /), rep(cl(iS u /))) to E; 
else add (rep(iS , u '),rep'(cl(iS' u '))) to E; 
for for each pair {S u >, S v >} in the MWSPD do 

if rep(5„/) does not have color c then add (rep(Sv), vep(S v ')) to E; 
else add (rep(<S' u '),rep'(<S'„/)) to E; 
end 

for each c-child u" of u' do 

if rep(cl(iSu/)) does not have color c then add (rep(5 u »), rep(cl(5 u '))) to E; 
else add (rep(5 u »), rep'(cl(5 u '))) to E; 
end 
end 



end 

26 end 

27 for each {S u , S v } in the MWSPD for which both u and v are multichromatic do 

28 if iep(S u ) and rep(S'^) have distinct colors then add (rep(5 n ), rep(5^)) to E; 

29 else add (rep(5 u ), Tep'(S v )) to E; 

30 end 

31 return the graph G = (S, E) 
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of nodes in the split-tree, which is 0(|»S|). 



□ 



Lemma 4.2 Let G be the graph computed by Algorithm [7} Let p and q be two points of S with 
different colors, and let {S u , S v } be the pair in the MWSPD for which p £ S u and q £ S v . Assume 
that u is a c-node for some color c. Then there is a path in G between p and rep{S u ) whose length 
is at most t'\pq\, where 



t' = 4\/d(/id + l)(l + 4/s) 



log (Vd(l + 4/a) 



+ 1, 



and s is the separation constant of the WSPD. 

Proof: Let w be the c-leaf such that p £ S w , and let w 
c-nodes that are on the path in T from w to u. 



u be the sequence of 



Recall from Definition 3.1 that each set S Wi , < i < k, has a representative iep(S Wi ) (of color 
c) associated with it. Also, recall the definition of the sets cl(S Wi ), < i < k; see Definition 3.1 
If cl(S Wi ) is a c'-node for some color c', then this set has one representative rep(cl(5 u , i )) associated 
with it. Otherwise, cl(S Wi ) is multichromatic and this set has two representatives rep(cl(S Wi )) and 
rep' (cl (£*„,.)) of different colors associated with it. We may assume without loss of generality that, 
for all < i < k, the color of rep(cl(S Wi )) is not equal to c. 

Let II be the path 



rep(cl(S W0 )) 
rep(cl(S UJ1 )) 



rep(S W0 ) 
rep(S wl ) 



rep(d(5 Wfc )) 



vep(S Wk ) = rep(5 n ) 
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The first edge on this path, i.e., (p, iep(cl(S Wo ))), is added to the graph G in lines 19-10 of the 
algorithm. The edges (rep(cl(S Wi )),rep(S Wi )), < i < k, are added to G in lines 
the edges (rep(S , l0j _ 1 ), 



15 Finally, 



It follows that II is 



, rep(cl(5 u , i ))), 1 < i < k, are added to G in lines 
a path in G between p and rep(5' u ). We will show that the length of II is at most t'\pq\. 

Let i be an integer with < i < k. Recall the definition of 01(5^); see Definition 3.1 We 
consider all pairs {S x ,S y } in the MWSPD, where x is a c-node on the path in T from W{ to the 
root, and pick the pair for which dist(S , a; , S y ) is minimum. We denote the pair picked by (S Xi , S Vi ). 
Thus, Xi is a c-node on the path in T from Wi to the root, {S Xi , S Vi } is a pair in the MWSPD, and 
cl(S Wi ) = S Vi . We define 

£i = dist(S Xi ,S y J. 

Consider the first edge (p,vep(cl(S WQ ))) on the path II. Since p £ S WQ C S XQ and vep(cl(S m )) £ 



S yo , it follows from Lemma 2.4 that 



\p,iep(d(S W0 ))\ < (1 + 4/a) • dist(^ , S yo ) = (l + 4/s)£ . 
Let < i < k and consider the edge (rep(cl(S lUi )), rep(5 U) J) on II. Since rep(S'^ i ) £ S Wi C S Xi and 



vep(cl(S Wi )) £ S Vi , it follows from Lemma 2.4 that 
(1) |rep(clOS Wi )),rep(S^)| < (1 + 4/a) • dist(^„ S m ) 
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Let 1 < i < k and consider the edge (rep(S , Wi _ 1 ), rep(cl(5 w J)) on II. Since rep(S , l0i _ 1 ) G S Wi _ 1 C S^.. 
and rep(cl((S' TOi )) 6 5^., it follows from Lemma 2.4 that 

IreptSV^reptcl^))) < (1 + 4/s) • dist^, S y J = (1 + 4/s)4 

Thus, the length of the path II is at most 

k 

£2(l + 4/a)£. 

8=0 

Therefore, it is sufficient to prove that 

k 

'Yjli < 2Vd(fj,d +1)(1 + 4/ 's) 2 \pq\. 

i=0 

It follows from the definition of cl(S u ) = cl(S Wk ) that 

4 = dist(5a; fc , 5 yfc ) < dist(S u ,S v ). 
Since, by Lemma [2~4| dist(5 u , S v ) < (1 +4/s)\pq\, it follows that 

(2) 4<(l+4/s)H. 

Thus, it is sufficient to prove that 

k 

(3) J2 £i <2v / d(M+l)(l + 4/s)4- 

i=0 

If = 0, then ([3]) obviously holds. Assume from now on that k > 1. For each i with < i < k, 
we define 

O'i = L max (f3(S Wi )) , 

i.e., aj is the length of a longest side of the bounding box of S w .. 
Let < i < k. It follows from Lemma 2A that 

2 

Lma,x(P{S Xi )) < — 4 
S 

Since m; is in the subtree of Xi, we have L max (f3(S Wi )) < L max ((3(S Xi ))- Thus, we have 
2 

(4) a { < -£i for < i < k. 



Lemma 12.61 states that 

(5) di < -o-i+d for < i < k — d. 

Let < i < k — 1. Since u>i is a c-node, there is a node w\ such that {S Wi ,S w >} is a pair in the 
MWSPD. Then it follows from the definition of cl(S Wi ) that 

£i = dist(S Xi ,S yi ) < dist(S Wi ,S w/i ). 
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By applying Lemma |2.7| we obtain 

dist(S Wi ,S w /) < d ( s + 4: ) L max (/3(5 , 7r(Wi ))) 
< — ^2 L max (/3(S Wi+l )) 

- 2 Oi+l- 

Thus, we have 

(6) £i < ^ Oi+i for < i < k - 1. 

First assume that 1 < A; < fid. Let < % < A; — 1. By using ([6]), the fact that the sequence 
OO) oi, . . . j ajfe is non-decreasing, and Q, we obtain 

Vd(s + 4) Vd(s + 4) r-, 
£i < a i+l < a k < Vd(l + 4/s)4- 

Therefore, 

< fcVd(l + 4/s)4 + 4 < (A; + + 4/s)4 < (lid + 1)^(1 + 4/s)4, 

8=0 

which is less than the right-hand side in ([3]). 

It remains to consider the case when k > fid. Let i > and j > be integers such that 
i + 1 + jd < k. By applying ^ once, ([5]) j times, and Q once, we obtain 



^ v^(s + 4) ^ ^(5 + 4) /l^ 

2 + 2 V 2 



a m+id < \/d(l +4/s) 4fl+jd- 



For j = fi = [log(\/d(l + 4/s))] + 1, this implies that, for < i < k — 1 - fid, 
(7) £ < hu+i+nd. 

By re-arranging the terms in the summation in ([3]), we obtain 

k fid \_{k-h)/(jid+\)\ 

i=0 h=0 j=0 

Let j be such that < j < [(k — h)/(fid + 1)J . By applying ([7]) j times, we obtain 



r -k-h-j(nd+l) < \ 2 j ^-fr- 



it follows that 

L(fe-h)/(AMi+l)J oo 



i=0 j=0 
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Thus, we have 



/id 



i=0 h=0 

By applying the fact that the sequence ao, oi, . . . , ajt is non-decreasing, followed by Q, we 
obtain, for < i < k — 1 and 1 < j < k — i, 



< 



Vd(s + 4) ^Vd(s + 4) r- 
— ^ Oj+i < y — a i+j < Vd(l + 4/ 



s)i 



i+j- 



Obviously, the inequality £i < \fd(l + 4/s)£j + j also holds for j = 0. Thus, for i = k — h and j = h, 
we get 

ik-h < V^(l + 4/s)4 for < /i < /id. 
It follows that 

^ 4 < 2 ^ + 4/s)4 = 2Vd(/id + 1)(1 + 4/s)4, 
j=0 h=0 



completing the proof that (|3]) holds. 



□ 



Lemma 4.3 Assuming that the separation constant s of the WSPD is chosen sufficiently large, the 
graph G computed by Algorithm^ is a t-spanner of the complete k-partite graph Kc x ...c k > where 
t = 2t' + 1 + 4/s and t' is as in Lemma 4-2. 



Proof: We denote the graph Kc x ...c k by K. It suffices to show that for each edge (p, q) of K, the 
graph G contains a path between p and q of length at most t\pq\. We will prove this by induction 
on the lengths of the edges in K. 

Let p and q be two points of S with different colors, and let {S u , S v } be the pair in the MWSPD 
for which p G S u and q £ S v . 



The base case is when (p, q) is a shortest edge in K. Since s > 2, it follows from Lemma 2.4 that 
u is a c-node and v is a c'-node, for some colors c and d with c ^ d . In line 17 of Algorithm [lj the 
edge (rep^u), rep(iSy)) is added to G. By Lemma 2.4 the length of this edge is at most (1+4/ s)\pq\. 



The claim follows from two applications of Lemma 4.2 to get from p to rep(5' u ) and from rep(5„) 
to q. 

In the induction step, we distinguish four cases. 

Case 1: u is a c-node and v is a c'-node, for some colors c and d with d. 
This case is identical to the base case. 

Case 2: u is a c-node for some color c and v is a multichromatic node. 



In lines 17-18 the edge (rep(5 u ), rep(5,;)) or (rep(S u ), rep(5(,)) is added to G. We may assume 



without loss of generality that (vep(S u ), rep(S v )) is added. By Lemma 2.4 the length of this edge 
is at most (1 + 4/s)\pq\. 



By Lemma 4.2, there is a path in G between p and rep(5' u ) whose length is at most t'\pq\. 
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First assume that q and iep(S v ) have the same color. Let r be a point in S v that has a color 



different from g's color. Since s > 2, it follows from Lemma 2.4 that \qr\ < \pq\. Thus, by induction, 
there is a path in G between q and r whose length is at most t\qr\, which, by Lemma 2.4 is at 
most (2t/s)\pq\. By a similar argument, since |r,rep(iS„)| < \pq\, there is a path in G between r 
and iep(S v ) whose length is at most (2t/s)\pq\. Thus, G contains a path between q and rep(S v ) of 
length at most (4t/s)\pq\. If q and rep(5„) have different colors, then, by induction, there is a path 
in G between q and rep(S v ) whose length is at most (2t/s)\pq\ < (4t/s)\pq\. 

Thus, the graph G contains a path between q and vep(S v ) of length at most (At/s)\pq\. 

We have shown that there is a path in G between p and q whose length is at most 

(8) (V + (1 + 4/s) + 4i/s) \pq\. 

By choosing s sufficiently large, this quantity is at most t\pq\. 

Case 3: u is a multichromatic node and v is a c-node for some color c. 
This case is symmetric to Case 2. 

Case 4: Both u and v are multichromatic nodes. 



In lines 28-29 the edge (rep(5 u ), rep(5„)) or (rep(5 , u ),rep(5(,)) is added to G. We may assume 



without loss of generality that (rep(5 u ), rep(S v )) is added. By Lemma 2.4 the length of this edge 
is at most (1 +A/s)\pq\. 

As in Case 2, the graph G contains a path between p and rep(5 n ) of length at most (4t/s)\pq\, 
and a path between q and rep(S*„) of length at most (4t/s)\pq\. 

It follows that there is a path in G between p and q whose length is at most 

(9) ((l+4/s) + 8t/s)\pq\. 

By choosing s sufficiently large, this quantity is at most t\pq\. □ 



Lemma 4.4 The running time of Algorithm^ is O(ralogra), where n= \S\. 

Proof: Using the results of Callahan and Kosaraju [3J, the split-tree T and the WSPD can be 
computed in 0(n log n) time. The representatives of all sets S u and all sets cl^u) can be computed 
in 0{n) time by traversing the split-tree in post-order and pre-order, respectively. The time for 



the rest of the algorithm, i.e., lines [3 ^31 is proportional to the sum of the size of T, the number of 
pairs in the WSPD and the number of edges in the graph G. Thus, the rest of the algorithm takes 
0(n) time. □ 



To summarize, we have shown the following: For any complete /c-partite geometric graph K 
whose vertex set has size n, Algorithm [T] computes a t-spanner of K having 0(n) edges, where t is 
given in Lemma 4.3 The running time of this algorithm is 0{n log n). By choosing the separation 



constant s sufficiently large, the stretch factor t converges to 
"1 



8v^ d 



logd 



+ d+l +1. 



In the next section, we show how to modify the algorithm so that the bound in Lemma 4.2 
is reduced, thus improving the stretch factor. The price to pay is in the number of edges in G, 
however, it is still 0{n). 
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5 An Improved Algorithm 



As before, we are given a set S of n points in M. d which is partitioned into k subsets C\, C%, . . . , Cf.. 
Intuitively, the way to improve the bound of Lemma [4 .2 1 is by adding shortcuts along the path from 
each c-leaf to the c-root above it. More precisely, from Q in the proof of Lemma 4.2 we know that 



if we go 1 + fid levels up in the split-tree T, then the length of the edge along the path doubles. 
Thus, for each c-node in T, we will add edges to all 25(1 + fid) c-nodes above it in T. Here, 5 is an 
integer constant that is chosen such that the best result is obtained in the improved bound. 

Definition 5.1 Let c S {1, 2, . . . , k}, and let u and v! be c-nodes in the split-tree T such that u' is 
in the subtree rooted at u. For any integer £ > 1, we say that u is C levels above v! , if there are 
exactly £ — 1 c-nodes on the path strictly between u and v! . We say that u' is a (-level c-child of u 
if u is at most ( levels above u' . 



The improved algorithm is given as Algorithm [2] The following lemma generalizes Lemma 4.2 



Algorithm 2: Computing a sparse (5 + e)-spanner of Kc\...c k - 



Input: A set S of points in R , which is partitioned into k subsets C\, . . . , and a real 
constant < e < 1. 

Output: A (5 + e)-spanner G = (S,E) of the complete fc-partite graph Kc x ...c h - 

Choose a separation constant s such that s > 12/e and (1 + 4/s) 2 < 1 + e/36 and choose an 
integer constant 5 such that < 1 + e/36. 

The rest of the algorithm is the same as Algorithm [T] except for lines 20-23 which are 
replaced by the following: 

for each ("-level c-child u" of v! do 



if rep(cl(5 u ')) does not have color c then add (rep(S , n "), rep(cl(5 u '))) to E; 
else add (rep(5 u »),rep'(cl(5 u '))) to E; 

if rep(cl(5 u «)) does not have color c then add (rep(cl(5 u //)), rep(5 u /)) to E; 
else add (rep'(cl(S u »)), rep(5 u /)) to E; 
end 



Lemma 5.2 Let G be the graph computed by Algorithm^ Let p and q be two points of S with 
different colors, and let {S u , S v } be the pair in the MWSPD for which p £ S u and q £ S v . Assume 
that u is a c-node for some color c. Then there is a path in G between p and rep(S u ) whose length 
is at most (2 + e/3>)\pq\. 

Proof: Let w be the c-leaf such that r G S w , and let w = Wq, w\, . . . , W}~ = u be the sequence of 
c-nodes that are on the path in T from w to u. As in the proof of Lemma |4.2| we assume without 
loss of generality that, for all < i < k, the color of rep(cl(S , u , 4 )) is not equal to c. 

Throughout the proof, we will use the variables Xj, yi, £j, and ctj, for < i < k, that were 
introduced in the proof of Lemma |4.2| 

We first assume that < k < 2S(fj,d + 1). Let IT be the path 

p -> rep(cl(S tl) )) rep(S u ). 
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It follows from Algorithm [2] that II is a path in G. Since p G S w = S Wo C S XQ and rep(cl(<S' w )) 



rep(cl(S' u , )) 6 S yo , it follows from Lemma 2.4 that 



(10) |p,rep(cl(SW)| < (1 + 4/a) ■ dist(S XQ , S yo ) = (1 + A/s)t Q . 

Since {S u , S v } is one of the pairs that is considered in the definition of cl(S WQ ), we have dist(5a; , S yo ) < 
dist(5 u , S v ). Again by Lemma [2~H we have dist(S u , S v ) < (1 + 4/s)\pq\. Thus, we have shown that 



|p,rep(cl(5 U) ))| < (l + 4/s) 2 H. 
By the triangle inequality, we have 

|rep(cl(5 11) )),rep(5' u )| < \rep(cl(S w )) , p\ + |p,rep(5 u ) 



Since p and rep(5 u ) are both contained in S u , it follows from Lemma 2.4 that \p,rep(S u )\ < 
(2/s)\pq\. Thus, we have 

|rep(cl( 1 S w )) ) rep(5 u )| < (1 + 4/s) 2 \pq\ + (2/s)\pq\. 

We have shown that the length of the path IT is at most 

(2(l+4/s) 2 +2/ S ) \pq\, 

which is at most (2 + e/3)|pg| by our choice of s in Algorithm [2} 

In the rest of the proof, we assume that k > 2S(p,d + 1). We define 

m = k mod (5(p,d + 1)) 

and 

. k — m 

777 = 

s{fj,d+iy 

We consider the sequence of c-nodes 

W = W , W6(nd+l)+mi w 28(nd+l)+m, ^3<5(^+l)+ mi ■■ ■ ,Wk = U, 

and define LT to be the path 

p -» rep(cl(5 W0 )) -> rep(5„,, (fid+1)+m ) 
-» rep(cl(5 W25(Md+1)+m )) -» rep(5 WM(Md+1)+m ) 
— ► rep(cl(S W34(M(i+1)+m )) — > rep(5 , W)34(M(J+1)+m ) 

-> rep(cl(5„,J) -> rep(S Wfc ) = rep(S u ). 

It follows from Algorithm [2] that IT is a path in C We will show that the length of this path is at 
most (2 + e/3)|pg|. 

We have shown already (see ([To|) ) that the length of the first edge on IT satisfies 



|p,rep(cl(5^))| < (l + 4/s)4. 
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2.4 



that 



The length of the second edge satisfies 

|rep(cl(^ )),rep(5^ (Md+1)+m )| < \rep(cl{S Wo )),p\ + \p,vep{S Wsifid+1)+ J\ 

< (l + 4/ S )4 + b,rep(5^ (/id+1)+m )|. 

Since p and rep(5^ i( ^ d+1)+m ) are both contained in S u , it follows from Lemma 

\p,vep(S Wd ^ d+1)+ J\ < (2/s)\pq\. 
Thus, the length of the second edge on II satisfies 

|rep(cl(5 u , )),rep(5 Si ,, (Md+1)+m )| < (1 +A/s)l + (2/s)\pq\. 
Let 2 < j < m! . We have seen in ([I]) in the proof of Lemma |4.2| that the length of the edge 

(rep(cl(5^ 5(Md+1)+m )),rep(S'^ a{M(i+1)+ra )) 
satisfies 

|rep(cl(S , WiS(Md+1)+m )),rep(^ j ., (M(i+1)+m )| < (1 + 4/s)£ j5{tld+1)+m . 
Again, let 2 < j < m! . Since 

rep(S UJ y_ :L ^ i5 ^ d+1 j +m ) £ S WjS ^ id+1 ^ +m C S x . s ^ d+ ^ +m 

and 

icep(cl(S WjS(iid+1)+m )) G Sy. S(jid+1)+m 
it follows from Lemma |2.4| that the length of the edge 
(rep(S" W(j _ 1) 

))) 

satisfies 

|rep(S^_ 1) , (Md+1)+m ),rep(cl(^ {M+1)+m ))| < (1 + 4/s)£ j5 ^ d+1)+m - 
We have shown that the length of II is at most 

( m' 
j6([id+l)+r 
3=2 

The definition of £q,£i, ■ ■ . ,£k implies that this sequence is non-decreasing. Thus, Iq < £s(fid+i)+m 
and it follows that the length of II is at most 

m' 

(2/ S )M + 2(l + 4/ S )^ 
i=i 

Recall inequality ([7]) in the proof of Lemma 4.2 which states that 
1 



h — 2^+^+1- 
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By applying this inequality 5 times, we obtain 

ii < (^) £i+5( M d+l)- 
For i = j5(fid + 1) + m, this becomes 



By repeatedly applying this inequality, we obtain, for h > j, 



(h-j)5 



l jS{^d+l)+m < ( 2 ) 4<5(/jd+l)+m- 



For h = m', the latter inequality becomes 

/] \ ( m '-j) s 
tjSivd+V+m < ( ^ ) 4- 



It follows that 

(m'-j)8 



in m f \ \ 

Zj6(nd+l)+m < Yl ( O ) 
rn'-l , 

til 



3=1 j=l 

rn'-l / -\ \ iS 

4 

v ^ / 

oo 

'•' ' \2> 
2 5 

= — 4- 

2 5 - 1 

According to ^ in the proof of Lemma 4.2 we have 

4 < (i + 4/ s )H. 

We have shown that the length of the path II is at most 
2/ S + 2(l+4/ S ) 2 2 ^M \pq\. 



Our choices of s and 5 (see Algorithm |2| imply that 2/s < e/6, (1 + 4/s) < l + e/36 and 
oiTzt < 1 + e/36. Therefore, the length of II is at most 

(e/6 + 2(1 + e/36) 2 ) \pq\ < (2 + e/3)\pq\, 

where the latter inequality follows from our assumption that < e < 1. This completes the proof. 
□ 
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Lemma 5.3 Let n = \S\. The graph G computed by Algorithm is a (5 + e)-spanner of the 
complete k-partite graph Kc\...C k an d the number of edges of this graph is 0{n). The running time 
of Algorithm mis 0(n log n) . 



Proof: The proof for the upper bound on the stretch factor is similar to the one of Lemma 4.3 



The difference is that instead of the value t' that was used in the proof of Lemma 4.3 
the value t! = 2 + e/3 of Lemma 



we now use 



5.2 



Thus, the stretch factor for the base case of the induction 



and for Case 1 is at most 

(1 + 4/s) + 2t' = 5 + 4/s + 2e/3, 

which is at most 5 + e, because of our choice for s in Algorithm [2] For Cases 2 and 3, the stretch 
factor is at most (see ^ in the proof of Lemma |4.3| where t = 5 + e) 

t' + (1 + 4/s) + 4i/s = 3 + e/3 + (4/s)(6 + e), 

which is at most 5 + e, again because of our choice for s. Finally, the stretch factor for Case 4 is at 
most (see ^ in the proof of Lemma 4.3 where t = 5 + e) 

(1 + 4/s) + 8t/s = 1 + (4/s)(ll + 2e), 

which is at most 5 + e, because of our choice for s. 

The analysis for the number of edges is the same as in Lemma |4.1[ except that the number of 
edges that are added to each c-node in the modified for-loop is 25{^d + 1) instead of one as is in 
Algorithm [II Finally, the analysis of the running time is the same as in Lemma 4.4 □ 



We have proved the following result. 

Theorem 5.4 Let k>2 be an integer, let S be a set of n points in M. d which is partitioned into k 
subsets Ci, C2, . . . , Cfc, and let < e < 1 be a real constant. In 0{n log n) time, we can compute a 
(5 + e)-spanner of the complete k-partite graph Kc 1 ...c k having 0{n) edges. 

6 Improving the Stretch Factor 

We have shown how to compute a (5 + e)-spanner with 0{n) edges of any complete /c-partite graph. 
In this section, we show that if we are willing to use 0{n log n) edges, the stretch factor can be 
reduced to 3 + e. We start by showing that a stretch factor less than 3, while using 0{n log n) 
edges, is not possible. 

Theorem 6.1 Let c > be a constant and let n and k be positive integers with 2 < k < n — 
2cy/nlogn. For every real number < e < 1, there exists a complete k-partite geometric graph K 
with n vertices such that the following is true: If G is any subgraph of K with at most c 2 n log n 
edges, then the stretch factor of G is at least 3 — e. 

Proof: Let D\, D2, and D3 be three disks of radius e/12 centered at the points (0,0), (l + e/6,0), 
and (2 + e/3,0), respectively. We place (n — k + l)/2 red points inside D\ and (n — k + l)/2 blue 
points inside D2. The remaining k — 2 points are placed inside D3 and each of these points has a 
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distinct color (which is not red or blue). Let K be the complete /c-partite geometric graph defined 
by these n points. We claim that K satisfies the claim in the theorem. 

Let G be an arbitrary subgraph of K and assume that G contains at most c 2 nlogn edges. We 
will show that the stretch factor of G is at least 3 — e. 

Assume that G contains all red-blue edges. Then the number of edges in G is at least 
((n — k + 2)/2) 2 . Since k < n — 2cy / nTogn, this quantity is larger than c 2 n log n. Thus, there 
is a red point r and a blue point b, such that (r, b) is not an edge in G. The length of a shortest 
path in G between r and b is at least 3. Since \rb\ < 1 + e/3, it follows that the stretch factor of G 
is at least , , 3 /Q , which is at least 3 — e. □ 



Theorem 6.2 Let k>2 be an integer, let S be a set of n points in M. d which is partitioned into k 
subsets Ci, C2, . . . , Cfcj and let < e < 1 be a real constant. In 0(n log n) time, we can compute a 
(3 + e)-spanner of the complete k-partite graph Kc 1 ...c k having O(nlogn) edges. 

Proof: Consider the following variant of the WSPD. For every pair {X, Y} in the standard 
WSPD, where \X\ < \Y\, we replace this pair by the \X\ pairs where x ranges over all 

points of X. Thus, in this new WSPD, each pair contains at least one singleton set. Callahan and 
Kosaraju jl] showed that this new WSPD consists of O(nlogn) pairs. 

We run Algorithm [2] on the set S, using this new WSPD. Let G be the graph that is computed 



by this algorithm. Observe that Lemma 5.2 still holds for G. In the proof of Lemma 5.3 of the 



upper bound on the stretch factor of G, we have to apply Lemma 5.2 only once. Therefore, the 
stretch factor of G is at most 3 + e. □ 



7 Conclusion 

We have shown that for every complete fc-partite geometric graph K in M. d with n vertices and for 
every constant e > 0, 

1. a (5 + e)-spanner of K having 0{n) edges can be computed in O(nlogn) time, 

2. a (3 + e)-spanner of K having O(nlogn) edges can be computed in O(nlogn) time. 

The latter result is optimal for 2 < k < n — 0(^/n logn), because a spanner of K having stretch 
factor less than 3 and having 0{n logn) edges does not exist for all complete A;-partite geometric 
graphs. 

We leave open the problem of determining the best stretch factor that can be obtained by using 
0{n) edges. 

Future work may include verifying other properties that are known for the general geometric 
spanner problem. For example, is there a spanner of a complete ^-partite geometric graph that has 
bounded degree? Is there a spanner of a complete /c-partite geometric graph that is planar? From 
a more general point of view, it seems that little is known about geometric spanners of graphs other 
than the complete graph. The unit disk graph received great attention, but there are a large family 
of other graphs that also deserve attention. 



18 



References 



[1] I. Althofer, G. Das, D. P. Dobkin, D. Joseph, and J. Soares. On sparse spanners of weighted 
graphs. Discrete & Computational Geometry, 9:81-100, 1993. 

[2] Prosenjit Bose, Paz Carmi, Mathieu Couture, Anil Maheshwari, Michiel Smid, and Norbert 
Zeh. Geometric spanners with small chromatic number. In Proceedings of the 5th Workshop 
on Approximation and Online Algorithms, Lecture Notes in Computer Science, Berlin, 2007. 
Spr inger- Verlag . 

[3] P. B. Callahan and S. R. Kosaraju. Faster algorithms for some geometric graph problems in 
higher dimensions. In Proceedings of the 4th ACM-SIAM Symposium on Discrete Algorithms, 
pages 291-300, 1993. 

[4] Paul B. Callahan and S. Rao Kosaraju. A decomposition of multidimensional point sets with 
applications to /c-nearest-neighbors and n-body potential fields. J. ACM, 42(l):67-90, 1995. 

[5] J. Gudmundsson and M. Smid. On spanners of geometric graphs. In Proceedings of the 10th 
Scandinavian Workshop on Algorithm Theory, volume 4059 of Lecture Notes in Computer Sci- 
ence, pages 388-399, Berlin, 2006. Springer- Verlag. 

[6] Giri Narasimhan and Michiel Smid. Geometric Spanner Networks. Cambridge University Press, 
New York, NY, USA, 2007. 

[7] Bhaskaran Raman and Kameswari Chebrolu. Design and evaluation of a new MAC protocol 
for long-distance 802.11 mesh networks. In MobiCom '05: Proceedings of the 11th annual 
international conference on Mobile computing and networking, pages 156-169, New York, NY, 
USA, 2005. ACM Press. 

[8] J. S. Salowe. Constructing multidimensional spanner graphs. International Journal of Compu- 
tational Geometry & Applications, 1:99-107, 1991. 

[9] P. M. Vaidya. A sparse graph almost as good as the complete graph on points in K dimensions. 
Discrete & Computational Geometry, 6:369-381, 1991. 



19 



